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P ! Abstract 

U ■ 

c/3 . We analyse the formation of cosmic structures in models where the dark matter 

^ ' is dominated by light gravitinos with mass of 100 eV - 1 keV, as predicted by gauge- 

mediated supersymmetry (SUSY) breaking models. After evaluating the number 
^ . of degrees of freedom at the gravitinos decoupling (gi^,), we compute the transfer 

^ I function for matter fluctuations and show that gravitinos behave like warm dark 

matter (WDM) with free-streaming scale comparable to the galaxy mass scale. We 
consider different low-density variants of the WDM model, both with and without 
cosmological constant, and compare the predictions on the abundances of neutral 
hydrogen within high-redshift damped Ly~a systems and on the number density of 
local galaxy clusters with the corresponding observational constraints. We find that 
none of the models satisfies both constraints at the same time, unless a rather small 
r^o value (^ 0.4) and a rather large Hubble parameter (^ 0.9) is assumed. Further- 
more, in a model with warm + hot dark matter, with hot component provided by 
massive neutrinos, the strong suppression of fluctuation on scales of ~ 1 /i~^Mpc 
precludes the formation of high-redshift objects, when the \ow-z cluster abundance 
is required. We conclude that all different variants of a light gravitino DM domi- 
nated model show strong difficulties for what concerns cosmic structure formation. 



^On leave of absence from Department of Physics, Tohoku University, Sendai 980-77, Japan. 



This gives a severe cosmological constraint on the gauge-mediated SUSY breaking 
scheme. 
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1 Introduction 



Since the moment (early 80's) that low-energy supersymmetry (SUSY) was invoked in 
gauge unified schemes to tackle the gauge hierarchy problem [|I|, it became apparent 
that it had also a major impact on several cosmological issues. By far the most studied 
consequence was the presence of a stable SUSY particle in all models where a discrete 
symmetry, known as R-parity, is imposed to prevent the occurrence of baryon and lepton 
renormalizable terms in the superpotential. Indeed, R-parity assigns a different quantum 
number to ordinary particles and their SUSY partners. Hence the lightest SUSY particle 
(LSP) is absolutely stable and constitutes, together with photons and neutrinos, a viable 
candidate for relic particles of the early Universe. 

The two best candidates we have to play the role of LSP are: the lightest neutralino 
(i.e. the lightest among the fermionic partners of the neutral gauge and Higgs fields) and 
the gravitino (the fermionic partner of the graviton in the gravity multiplet) 0. Which 
of the two is the actual LSP strictly depends on the mechanism one envisages for the 
SUSY breaking, or, more precisely, for the transmission of the breaking of SUSY from 
some hidden sector to the observable sector of the theory (ordinary particles and their 
superpartners belong to this latter sector). If the "messengers" of the SUSY breaking are 
of gravitational nature (as it happens in the more "orthodox" supergravity models), then 
the lightest neutralino is likely to be the LSP. In these schemes the gravitino mass sets 
the scale of SUSY breaking in the observable sector and, hence, it is expected to be in 
the 10^ — 10'^ GeV range. On the other hand, it has been vigorously emphasized recently 
(after ten years of silence about this alternative) that gauge, instead of gravitational 
interactions may be the vehicle for the transmission of the SUSY breaking information 
to the observable sector . In these scenarios the scale of SUSY breaking is much lower 
than in the supergravity case and consequently, as we will see below, the gravitino mass is 
much lower than 10^ GeV. Hence in this class of gauge mediated SUSY breaking (GMSB) 
models the gravitino is more likely to play the role of LSP with a mass which can range 
a lot, depending on the specific scale of SUSY breaking, say from a fraction of eV up to 
O(GeV). 

/^From a cosmological point of view, the neutralino LSP scenario with a lightest neu- 
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tralino in the tens of GeV range constitutes an ideal ground for a cold dark matter (CDM) 
proposal 0]. Indeed there exists a sufficiently vast area of the SUSY parameter space 
where such an LSP becomes non relativistic at a sufficiently early epochs so as to make 
its free-streaming mass much smaller than the typical galaxy mass scale (~ lO^^AfQ). The 
standard version of the CDM scenario, with Qq = 1 for the density parameter, h = 0.5 
for the Hubble parameter^ and P{k) oc k for the post-inflationary power spectrum of 
Gaussian adiabatic density fluctuations, is generally accepted to fail in reproducing sev- 
eral observational tests. On scales of few tens of h~^Mpc it develops a wrong shape of the 
power spectrum ||^. Furthermore, once normalized to match the detected level of CMB 
temperature anisotropics 0, it produces too large fluctuations on scales ^ 10 /i~^Mpc, 
with a subsequent overproduction of galaxy clusters [0. 

These failure of the standard CDM model may be overcome in these SUSY models 
by flnding a way to suppress fluctuations on 10 h^^Mpc scales, without decreasing too 
much power on the ~ 1 h^^Mpc scale, which would delay too much the galaxy formation 
epoch. A flrst possibility is adding to the LSP CDM candidate some massive light neutrino 
(Cold+Hot DM model; CHDM) to provide about 20-30% of the critical density §. This 
has just the effect of decreasing the fluctuation amplitude around the neutrino free- 
streaming scale, so as to change the power-spectrum in the right direction. A further 
possibility is assuming a density parameter substantially smaller than unity, either with 
or without a cosmological constant to provide spatial flatness . A lower cosmic density 
gives rise to a larger horizon size at the matter-radiation equality epoch, so as to increase 
the large-to-small scale power ratio in the spectrum of cosmic density fluctuations. 

If the gravitino is the LSP, one loses the traditional CDM candidate, being such 
a gravitino a more likely warm (WDM) candidate, its free-streaming mass scale being 
comparable to the galaxy mass scale 0, |ll] . This happens when its mass lies in the range 
[0.1-1] keV, which represents the situation that we will analyse in detail in this paper. 

It is already known that just replacing the cold LSP with a warm one in the standard 
CDM scenario does not provide a viable scenario for the formation of cosmic structures 
||12|| . Indeed, the effect of introducing the warm component is that of suppressing fluctu- 
ations only at the galaxy mass scale, while leaving the power spectrum unaffected on the 
^We take Hq — Mpc^^ for the Hubble constant. 
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cluster mass scales, where standard CDM fails. 

Therefore, if we desire a GMSB scheme to provide the dominant DM content of the 
Universe, we need some prescription to improve the WDM scenario. To this purpose, we 
will analyse in the following what happens if we follow the same pattern as for improving 
CDM, namely either adding a hot neutrino component or lowering the density parameter. 
Our analysis will focus on the interesting class of GMSB schemes, although many of our 
conclusions may equally well apply to models with a generic WDM other than the light 
gravitino. 

The purpose of our analysis is twofold. On one hand, given the success of suitable 
CHDM and low-density CDM models in accounting for several observational constraints 
(in particular providing a low level of density fluctuations at the 10/i~lMpc scale to avoid 
cluster overproduction, while having enough power at about l/i~^Mpc to form galaxies at 
an early enough epoch) , we ask whether the agreement can be kept when a warm gravitino 
component replaces the cold candidate. On the other hand, from a more particle physics 
oriented point of view, we would like to make use of the cosmological constraints related 
to the DM issue to infer constraints on the GMSB models, in particular shedding some 
light on the range of the allowed (or at least cosmologically favoured) scales of SUSY 
breaking in this class of theories. 

The paper is organized as follows. In Section 2 we present the general features of the 
GMSB models, focusing in particular on their predicted light gravitinos. We compare the 
two scenarios, gravity-mediated and gauge-mediated SUSY breaking, in relation to their 
LSP predictions and imphcations for DM. We provide the main tools for the computation 
of the rehc gravitino abundance in GMSB models. Section 3 describes the scenarios for 
the formation of cosmic structures when the DM content is dominated by light gravitinos. 
Here we compute the corresponding power spectra of density fluctuations at the outset 
of recombination. Afterwards, we present the observational data that we will use to 
constrain this class of models, namely the abundance of neutral hydrogen within high- 
redshift damped Ly-a systems and the number density of local galaxy clusters. In Section 
4 we compare the model predictions for the formation of cosmic structures with the 
abovementioned data. The main conclusions of our analysis are summarized in the final 
Section 5. 
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2 Light gravitinos in SUSY 

In a supersymmetric model each ordinary particle is associated with a superpartner. 
We assign R-parity even to the ordinary particles and odd to their superpartners. In 
supergravity, that is a natural extension of the supersymmetric standard models to the 
framework of local supersymmetry, we have another R-odd particle, the gravitino, which 
is the superpartner of the graviton. The lightest of the R-odd particles, namely the 
lightest superparticle (LSP), is absolutely stable, under the assumption of the R-parity 
conservation, which was originally introduced in order to avoid too fast proton decays. 
The LSP is thus a dark matter candidate, if its expected relic abundances lie within a 
suitable range of values. 

As a starting point, we review some properties of the gravitino. Imposing the vanishing 
cosmological constant in the Einstein supergravity Lagrangian, one finds that the gravitino 
mass is related to the SUSY breaking scale Asusy as follows: 

(1) 

where Mpi is the reduced Planck mass ~ 2.4 x 10^^ GeV. On the other hand, the soft 
SUSY breaking masses for the superparticles are given as 

^soft J^, [2) 

where M effectively represents the mass scale of the interactions that transmit the break- 
down of SUSY in the hidden sector to the observable sector, the latter including particles 
of the SUSY standard model. We call M the messenger mass scale. In the conventional 
scenario of the gravity-mediated SUSY breaking, the transmission is due to gravitational 
interactions. In this case, the messenger mass scale is M ~ Mpi, so that the gravitino 
mass will be comparable to the other soft masses. In order to have the soft masses 
at the electro-weak scale the SUSY breaking scale should be at an intermediate scale 

On the other hand, one can consider the case where the SUSY breaking is transmitted 



by gauge interaction. The idea of the gauge mediation |]T3[ is older than the gravity- 



mediation, and has recently been revived with fruitful results [H]. In this case the gauge 
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interaction can set the messenger mass scale much lower than the Planck mass. Since 
the soft masses are fixed at the electro-weak scale, the SUSY breaking scale can be much 
smaller than the intermediate scale of ^Jm^Mpi. Correspondingly the gravitino can be 
much lighter than the other superparticles. Now a crucial question is: how light is the 
gravitino? The answer should depend on the details of the messenger of the SUSY break- 
ing. In most of the gauge-mediated models, there are three independent sectors. They 
are the hidden sector, the messenger sector and the observable sector. The interaction 
between the last two sectors is the standard-model gauge interaction, so its strength is 
fixed. But the interaction between the first two is model dependent, so is the messenger 
mass scale. For example, in the original models of gauge-mediation it was shown 



that h^susY cannot be smaller than 10'' GeV, the corresponding gravitino mass being 
~ 10^ keV. However, a lighter gravitino should be possible from viewpoints of both model 
building and phenomenology. In the SUSY gauge-mediated approach the soft masses arise 
at the loop level (to avoid the supertrace constraint [1^). Hence a lower bound on Ksusy 
is provided by the relation 

^ 167r^ 

J^SUSY^ ^2 ^soft (3) 



where g is some gauge coupling constant. For instance, recently Izawa et al. [0 have 
constructed a model where rrisoft^ 0.1(7^/167r^A5[/5y. In this case, Asusy can be as small 
as O(IO^) GeV for nisoft = O(IO^) GeV. In view of the above consideration, in this paper 
we consider the following gravitino mass range Q 

1 eV < m^< a few TeV (4) 

It is noteworthy that interaction of the longitudinal component (spin 1/2 component) 
of the gravitino is fixed by the low-energy theorem. Namely the would-be Goldstino has a 
derivative coupling to the supercurrent with l/Ag^^y = I/v^^^g^p; suppression. After 
integration by parts and the use of equations of motion 0, we obtain the following effective 
■^In the framework of no-scale models, one may consider a somewhat larger range of the gravitino 



masses |17 . 

■*Here we present the formulas for massless gauge bosons. One needs to modify the formulas when the 
gauge bosons get massive due to symmetry breakdown. In our numerical computation in Section 2.2 this 
correction is taken into account. 
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Lagrangian |T^, 

2 2 

>Ce// = o ^[7.. 7.] AF,. + ^^^dxL^ + /^.c, (5) 

where G represents the longitudinal component of the gravitino (the Goldstino) and m^, 
and are the masses of a gaugino A, a chiral fermion x ^i-nd its superpartner 0, 
respectively. The point is that as the gravitino mass gets smaller the interaction becomes 
stronger. What happens physically is that a lighter gravitino corresponds to a lighter 
messenger scale, and therefore the Goldstino which is in the hidden sector has a stronger 
interaction to the fields in the observable sector. This point is crucial when we discuss 
the cosmology of the light gravitino. 

2.1 Two Scenarios: Neutralino LSP and Gravitino LSP 

Among the superparticles which appear in the supersymmetric standard models, a neu- 
tralino tends to be the lightest one and, therefore, it is stable. The neutralino LSP with 
mass of the order of 100 GeV turns out to be a good candidate for the cold dark matter 
(CDM) p[. In gravity-mediated models with ~ 10^-10'^ GeV, we face the traditional 



gravitino cosmological problems |T9[. Namely, unless gravitinos are strongly diluted at 
inflation and they are not regenerated in the reheating phase {Treh^ 10*^ GeV), they would 
spoil the canonical picture of big-bang nucleosynthesis (BBN). 

On the other hand, if the gravitino is lighter than the neutralino, the latter is no longer 
stable, and decays to the gravitino. It was pointed out that its decays would also 
destroy the BBN if its life time is sufficiently large. A limit on the life time depends on 
the abundances of the neutralinos before decay. We quote here a conservative bound of 
10^ sec as an upper bound for the life time of the neutralino from the BBN constraint. 

In this case the gravitino will be the stable LSP. Suppose that the spin 1/2 compo- 
nents of gravitinos were in thermal equilibrium at an early epoch.| As temperature went 
down, the processes which kept the gravitinos in equilibrium became ineffective and they 
decoupled from the thermal bath. After that, the number of gravitinos per comoving 
volume was frozen out. This freeze-out took place while the gravitinos were relativistic. 



^We assume that the Universe underwent inflationary era, so that the spin 3/2 components of the 
gravitinos were not thermahzed after that. 
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Following a standard procedure one can calculate the relic density of the gravitinos 



22 



/100\ 



(6) 



where Qq is the contribution of the (thermal) gravitinos to the density parameter, h is 
the Hubble parameter in units of 100 km/s/Mpc, and (?* stands for the effective degrees 
of freedom of relativistic particles when the freeze-out of the gravitinos takes place. Note 
that = 106.75 for the full set of particle contents of the minimal standard model 
and (?* = 228.75 for those of the minimal supersymmetric standard model. Thus one 
expects that (7* at the freeze-out will fall somewhere in between the two numbers. The 
computation of is a crucial point of our analysis and we will come back to it later on. 
For later convenience, we introduce the yield, Y^o, of the gravitinos, defined by 

\ s / 00 g* 

where is the number density of the gravitinos and s is the entropy density. The sub- 
script 00 means that the ratio is evaluated at a sufficiently late time (i.e., low temperature) 
at which it is constant. 

We will first briefiy discuss the case when the relic abundance of the gravitinos cal- 
culated in this way exceeds the closure limit, 1. This corresponds to the gravitino 
mass region m^K^ 1 keV ((7*/100)/i^. In this case, as was discussed in Refs. H . 



en- 



tropy production is needed to dilute the gravitino abundance in order not to overdose the 
Universe. To avoid an excessive reproduction of the gravitinos after the entropy produc- 
tion, its reheating temperature must be low; its upper bound varies from 10^ to 10® GeV, 
depending on the gravitino mass. The lower the gravitino mass is, the lower the reheat- 
ing temperature should be. If the reheating temperature happens to saturate the upper 
bound quoted above, the gravitinos will dominate the energy density of the Universe, and 
play the role of DM. 

On the other hand, the low reheating temperature required by the closure limit leads 
to the question of how to generate the baryon asymmetry of the Universe. Since the 
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reheating temperature can be still higher than the electro- weak scale, baryogenesis during 



the electro-weak phase transition may work for some region of the parameter space p3 |. 
Another possibility is to use the Affleck-Dine mechanism, which was explored in detail in 
Ref. in the framework of the gauge-mediated SUSY breaking. 



When the gravitino mass is smaller than {g^/100)h'^ keV, the thermal relic density 
of the gravitinos Qq is smaller than one. This is the region that we will study in detail 
in this paper. As we discussed previously, models providing this range for the gravitino 
mass can be devised. It is also interesting to mention that a possible explanation of the 
ee77 event at CDF by the light gravitino scenario requires this range of gravitino 



mass; otherwise the neutralino would not decay into a photon and a gravitino inside 
the detector. A particularly interesting parameter region for cosmology is the region in 
which 0.1^ 1 is realized, and thus the gravitino mass density constitutes a significant 
portion of the density of the whole Universe. A DM particle with mass within the sub- 
keV to keV range is known as warm dark matter [p6| , p!0| , |Tl[. Differently from CDM, it 
is characterized by having a sizable free streaming length until matter-radiation equality, 
roughly of the order of Mpc, but still much smaller than that of the hot dark matter 
(HDM), like a few eV neutrino. We will discuss scenarios of cosmic structure formation 
within a WDM dominated Universe in the following sections. 

If, instead, the gravitino mass is as small as to give Qq <^ 0.1, then it becomes 
cosmologically irrelevant and an alternative DM candidate is required. 

2.2 Computation of 

Before moving to the discussion of cosmic structure formation, we would like to come back 
to the question of the effective degree of freedom of relativistic particles at the freeze- 
out of gravitinos. Of particular interest is the region where 1 keV so that gravitinos 
of thermal origin dominate the energy density of the Universe. The crucial relevance of 
lies in the fact that, for a specified value of VLqU"^ , it fixes the corresponding gravitino 
mass and, therefore, the free-streaming scale. 

The production and destruction rates of the gravitinos due to scattering processes are 
proportional to the fifth power of the temperature and, therefore, their abundance rapidly 
drops down as the temperature decreases. Thus, decay and inverse decay processes are 
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more important for a light gravitino whose freeze-out occurs at a rather low temperature 



20| . In the following we will focus on these processes. 

The relevant Boltzmann equation can be casted in the form 

+ mn^ = C, (8) 

where Uq is the gravitino number density and H is the expansion rate of the Universe. 
As a collision term, we consider contributions from two body decay (and inverse decay) 
processes 

(9) 



C = Er(a^6G)(|^)n. (^1 - 



Here r(a hG) is the partial width of the species a to 6 and G, {nia/Ea) stands for the 
thermal average of the Lorentz boost factor, with ma and Ea being mass and energy of a, 
Ua is its number density and finally the superscript "eg" indicates the equilibrium value 
of a given quantity. After some algebra, the above Boltzmann equation can be rewritten 
as 



s' , s' 



Y' RY = RY^i, (10) 

Y.T{a^hG){ma/Ea)na/n^ 
^ = ) (-'--Lj 

where Y is the yield of the gravitinos as defined by Eq. (|^, and the apex symbol denotes 
derivative with respect to the temperature. Eq. (|T0[) can be solved to give 

Y{T) = Y^'^iT) + J^" dT'exp dT" R{T")s' /3s^ Y^" '{T'). (12) 

Here the temperature Tq is taken to be sufficiently high so that the gravitino is still in 
thermal equilibrium. 

In order to understand the meaning of Eq. (|12|), let us consider the case where R{T) 
changes abruptly at a temperature Tj such as R{T) = oo for T > Tf and for T < Tf. In 
this case we can approximate exp(— /J dT" R{T")s' /3s) with a step function 9(Tf — T'), 
so that 

Y{T) = r^«(T) + r° dT'eiTf - T')Y''^ '{T') = Y^\Tf), (13) 
Jt 
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thus reproducing the usual result. In the present case, however, R(T) gradually decreases 
as a species becomes non-relativistic. Therefore, we need to integrate Eq. (|12D numerically 
to evaluate Y{T) accurately. For sufficiently low T, the yield Y{T) approaches its constant 
value Yoo, from which we obtain using Eq. (|^). 

Results are presented in Table 1. We show the value g^: for a range of model parameters. 
In this computation, we assumed a typical sparticle mass spectrum in a simple class of 
gauge- mediated models p^. Explicitly, we take for the gauginos 



M, = l^Aa, Ak = ^Aa, M, = ^Aa, (14) 



and for the sfermion masses 



= 2 



^- (15) 



In the above expressions is a gauge coupling constant in the standard model, Y is 
a hypercharge of f/y(l), while C3 = 4/3 for a SU{3)c triplet, C2 = 3/4 for a SU{2)l 
doublet, and otherwise. A^, A5 are introduced to parameterize the transmission of 
SUSY breaking from the messenger sector to the observable sector .0 We provide (7* values 
for two cases: (a) the right-handed slepton mass m^-^ equals to the bino mass Mi, i.e. 
rrij^ = Ml, and (b) rrij^ = 2Mi. In both cases, we find that g^, is around 100 for a 
wide range of the parameter space. For a given Q^, a lower value of g^^ implies a lighter 
gravitino, making structure formations at small scales more difficult, as we will discuss in 
the following sections. The fact that tends to lie in the lower side should be kept in 
mind, though we will explore a somewhat wider range for g^. 

Table 1: Value of effective degrees of freedom of relativistic particles, at the gravitino 
freeze-out, function of the gravitino mass rriQ and the f/(l)y gaugino mass Mi. 

In the case (a) the right-handed slepton mass is mj^ = Mi, and in (b) mj^ = 2Mi. 

3 Light gravitinos and cosmic structure formation 



^To avoid further complication, we set a light Higgs mass to be the Z° mass, and masses of heavier 
Higgs and higgsinos to be the same as the left-handed slepton mass. Furthermore we did not include 
D- or F-term contributions to the scalar masses. Also we ignored the mixing in the mass matrix of the 
neutralino and the chargino sector. 
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3.1 Computation of the transfer functions 

The fundamental quantity that allows to make predictions about the formation of cos- 
mological structures, once the underlying Friedmann background is fixed, is the transfer 
function T{k), which convey all the informations about the evolution of a density fluctua- 
tion mode at the wavenumber k through the matter-radiation equality and recombinations 
epochs. In the following we will discuss how the transfer functions for the models under 
consideration are computed. As for models containing only the warm gravitinos (WDM) 
we will consider the fio ^ 1 cases, both with (ACDM) and without (OCDM) a cosmo- 
logical constant term, = 1 — Qq, to restore the spatial flatness. Furthermore, we will 
consider also the class oi Qq = ^ mixed models, whose DM content consists both of warm 
gravitinos and one species of hot neutrinos, having mass ^ 91 Qi^h"^ eV {Qi, is the 
neutrino contribution to the density parameter). 

Here we will only sketch our implementation of the Boltzmann code to compute T{k) 
and we refer to the relevant literature ([^]; |^) for more technical details. 

The transfer function is defined as 

Of)- 



T{k) 



■ ,z=0 



Of)- 



(16) 



where Ns is the number of different massive species in the model, 6i is the energy over- 
density of the i — th component and Zi a suitable initial redshift such that the smallest 
considered scale is much larger than the horizon scale at Zi. 
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We evaluate the transfer function for the models of interest in two steps: firstly we solve 
the equations for the fluctuation evolution of all the species involved in the models (namely 
the baryons, the radiation, the massless and massive neutrinos, and the gravitinos) for 
a number of k values; secondly, we find a suitable analytic expression which is able to 
provide a good fitting to the transfer functions for the whole class of considered models, 
by varying a minimal set of parameters. 

As for the fluctuation evolution, the goal is to find the final amplitude 6i for the 
different species, given the initial one. This goal is achieved in different ways for different 
components. For baryons only two differential equations must be solved: one regarding 
their overdensity and one for their velocity; for relativistic particles it is necessary to solve 
a hierarchy of coupled differential equations for the coefficients of the harmonic expansion 
of the perturbation in order to well describe the free-streaming behaviour . 

For massive free-streaming particles, different free-streaming behaviours can be ex- 
pected depending on which fraction of particles has to be considered still relativistic at 
a certain epoch. For these species it is therefore necessary to follow the fluctuation evo- 
lution separately for particles having different momenta. A representative set of different 
values of the momentum is chosen, and the density fluctuation evolution is evaluated for 
each value of this set. The overall 6i is therefore found by integrating the zero-th order 
harmonic coefficients over the momentum, with weights chosen on the basis of the dis- 
tribution function. This is the reason why, unlike for CDM, for massive free-streaming 
components the shape of the spectral distribution function affects the shape of the final 
transfer function. 

In our case, both gravitinos and massive neutrinos have an initial thermal distribu- 
tion, so the equations describing their evolution are qualitatively the same for both the 
components. What makes the difference between the two is the redshift at which they 
become non-relativistic, being higher for the warm G than for hot u. As a consequence, 
such two particle populations will be characterized by different free-streaming scales. 

All the calculations were performed in the syncronous gauge. For a detailed description 
on how a thermal free-streaming component is treated in the syncronous gauge, see Mah 
and Bertschinger [5S|. From a numerical point of view, we find that a higher degree of 
accuracy is needed when dealing with WDM-dominated models if compared to the CDM- 
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dominated ones. The reason is that all the Si are coupled by means of the potential; whose 
evolution equation, in turn, depends upon all the the different over densities, each of them 
contributing with a weight Qi. If the overdensity of the most abundant component is not 
well evaluated, the error propagates via the potential to all the other components, and 
over time. In the case of standard MDM, CDM plays this role, it stabilizes the value 
of the potential so that a lower accuracy in the integrals over the momenta of the hot 
component is allowed. 

In the models considered hereafter, gravitinos and massive neutrinos are the most 
abundant components, and their overdensities are evaluated by mean of integrals. It is 
therefore necessary to choose the integration method that, at the same time, (i) provides 
the best accuracy, and (ii) minimizes the number of values of the momentum over which 
the integration is performed, so as to keep the number of differential equations to be 
solved as small as possible. 



Within the class of Gauss integration methods we verified that, keeping fixed the 
number of integration points in momentum space, Gauss-Legendre integration performs 
better than Gauss-Laguerre, especially for high values of k . Furthermore, we found that 
using Gauss-Legendre integration, 20 integration points are adequate to obtain stable 
results. 

We computed the transfer function up to k^ax — 1 Mpc~^ (for fio = 1 and h = 0.5), 
with higher k values requiring too high an accuracy to be reached within a reasonable 
computational time. We will show in the following that such a kmax value is larger than the 
free-streaming wavenumber, kfg. Therefore, we expect that the behaviour of the transfer 
function at A; > k^ax has a marginal influence on the hierarchical clustering regime at 
k < kfs, we are interested in. 

In order to provide an analytical fitting to the transfer functions for the class of purely 



WDM models, we resorted to the expression provided by Bardeen et al. |^| 

TwDMik) = TcDMik) exp f-^ - ih^] , (17) 



where 



Tr 



CDM 



ln(l + 2Mq) 
2.34g 



1 + 3.89g + {IG.lqf + (5.46^)^ + (6.71g)- 



-1/4 
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is the transfer function for CDM models. Here, q = k/Th and the expression for the shape 
parameter, T = Q,Qhexp(—ilB — V2hilB/^o) accounts for the presence of a non-neghgible 
baryon fraction Qb P^ - 

Therefore, by fitting the transfer function, as computed by the Boltzmann code, with 
eqs.([T7|) and (|18D one obtains the value for the free-streaming scale, Rfg- More in detail, 
our procedure to estimate Rfs proceeds as follows. 

(a) We run the Boltzmann code assuming Qq = 1 and taking (7^ = 100 and 200; the first 
value is rather representative of realistic cases, while the larger corresponds to a 
very cold G population. 

(b) The free-streaming scale for the fio < 1 cases is then computed by resorting to the 
scaling relation Rf^ oc oc Qq (cf. eq.(^), where Qq = Qq — Qb- 

As a result, we find that 

/ a \-4/3 

Rfs = 0.51 {nchV [j^) Mpc (19) 

always provides an accurate fitting of the exponential suppression of fluctuations on small 
scales. We note that our value for Rfs is larger by a factor ~ 2.5 than that given by 



Kawasaki et al. This difference mainly comes from the fact that our value is directly 
obtained by fitting the exactly computed transfer function, while their value comes from 
the usual relation between Rfs and Znr (see, e.g., eq.(9.88) in the Kolb & Turner book [pi[]). 



the redshift at which gravitinos becomes non relativistic, that represents an approximation 
to the Rfs value. We also confirm the warning by Bardeen et al. who pointed out 
that the exponential cutoff in eq.([17D marginally underestimates the transfer function 
on intermediate scales, 0.1^ 0.5 {Qoh'^)^^ Mpc~^. However, we did not attempt here 
to look for a more accurate fitting expression, since (a) the effect is always quite small 
(^ 5-10%) and (b) we will mainly concentrate our analysis on the small scales relevant 
to galaxy and galaxy cluster formation. 

We plot in Figure 1 the TwDM^k) shape for i^o = 1 ^oi different values (left panel) 
and for two fio < 1 cases (right panel), also comparing with the corresponding CDM 
cases. It is apparent the power suppression on small scales, which depends both on g^: 
and on the parameters of the Friedmann background (cf. eq.(p!9|)). 
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As for the warm + hot DM (WHDM) case, transfer functions have been computed 
for Qi, = 0.1,0.2,0.3,0.4 and 0.5 in the case of only one massive neutrino (cf. ref. |^ 
for the effect of introducing more than one massive u), taking = 100 and 200 and 
always assuming fio = 1. The analytical fitting is provided by eg. (P!7|) , where the CDM 
transfer function is replaced by the CHDM one, as provided by Pogosyan & Starobinski 



33| . Taking = 1 — iliy — Hb, we find that eq.(0) always provides an accurate fitting 



to the exponential cutoff in the transfer function. The shapes of Ty/HDKiik) are plotted 
in Figure 2, showing both the effect of changing at fixed VLy (left panel) and the effect 
of changing Vt^, at fixed g^ (right panel). 

According to Figs. 1 and 2, it turns out that the effect of replacing the CDM compo- 
nent with light gravitinos of mass given by eq.(^) is that of eliminating the hierarchical 
clustering below some free-streaming mass scale. In order to provide an estimate of the 
free-streaming mass scale, we resort to the almost Gaussian cutoff at large fc, to define it 
as 



Mf, = {2TiR%Y"-p ^ 0.55 (^^j in^h')-'^oh'M,2, (20) 

where p is the average cosmic density and M12 = IO^^Mq. Therefore, eq.(pOD provides the 
limiting mass for the development of hierarchical clustering: structures of smaller masses 
form after structure of mass larger than Mfg, as a product of their fragmentation. As a 
consequence, we expect that a crucial constraint for the whole class of WDM-dominated 
models will come from the abundance of high-redshift cosmic structures. 

Having fixed the expression for the transfer function, we define the power spectrum 
of the density fluctuations as P{k) = AT'^{k)k^p'~ , where Upr is the primordial (post- 
inflationary) spectral index. The amplitude A is determined by following the recipe by 
Bunn & White to normalize both low-density fiat and open models to the 4-year 
COBE data. In the following, we will not consider the case of non-negligible contribution 
of tensor mode fluctuations to the CMB anisotropies. Indeed, such an effect would lead 
to a smaller spectrum amplitude, with a subsequent delay of the galaxy formation epoch 
that, as we will see, represents a major problem for WDM-dominated models. 

We plot in Figure 3 the r.m.s. mass fluctuation o"a/ for the same models whose T[k) 
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have been plotted in Fig. 1. This quantity is defined as 

crl = ^ I W\kRM) , (21) 

where the length scale associate to the mass scale M, Rm = {inp/Sy^M^/^, is the 
radius of the top-hat sphere whose Fourier representation is given by W{x) = 3(sinx — 
a;cosx)/a;^. For each model, the corresponding free-streaming mass scale corresponds 
to the transition from heavy to light curves in Fig. 3, while the completely light curves 
represent the corresponding CDM cases. It is apparent that such a scale is always at 
least of the order of a large galaxy halo. The flattening of cxm at small masses represents 
the imprint of non-hierarchical clustering. On the other hand, it turns out that the 
behaviour on the scales of galaxy clusters, ~ IO^^H'^Mq, is rather similar as for the 
CDM-dominated case. In the following we will use the abundance of local galaxy clusters 
and of high-redshift protogalaxies, through data about damped Ly-a systems, to constrain 
the whole class of WDM-dominated models. Constraints on larger scales, like bulk-flows 
data are much more similar to the CDM case. 

3.2 Observational constraints 
3.2.1 High-redshift objects 

The flrst constraint that we consider comes from the abundance of neutral hydrogen (HI) 
contained within damped Ly-a systems (DLAS; see ref. for a review about DLASs). 
DLAS are observed as wide absorption through in quasar spectra, due to a high HI column 
density (^ 10^° cm~^). Since at 3 the fractional density of neutral hydrogen associated 
with DLASs, Qhi, is comparable to that associated to visible matter in local galaxies, 
it has been argued that DLASs trace a population of collapsed protogalactic objects. 
In this context, a crucial question is to understand whether the observed Qhi provides 
a fair representation of the collapsed gas fraction at a given redshift. Effects like gas 
consumption into stars, ampliflcation biases due to gravitational lensing of background 
QSOs [^] and dust obscuration ||38| could well alter flnal results. However, such effects 



are believed to play a role at low redshift {z ~ 1-2), while they are expected to be less 
relevant at the highest redshifts at which DLAS data are available. For this reason, we will 
consider as the most constraining datum the value of Qhi reported by Storrie-Lombardi 
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et al. ||39| at redshift z ~ 4.25 and will assume that all the HI gas at that redshift is 
involved in the absorbers. 

Several authors recognized DLASs as a powerful test for DM models using both linear 
theory and numerical simulations |^D|. The recent availability of high-resolution spectra 
for several DLAS systems, allowed Prochaska & Wolfe ^l|] to use the internal kinematics 
of such systems to severely constrain a CDM model. 

In order to connect model predictions to observations, we consider the fraction of DM 
which at redshift z is collapsed into structures of mass M, 



Qcoii = erfc 



(22) 



Accordingly, Qhi = ^b^coU- Here, aM^z) is the r.m.s. fluctuation at the mass scale M 
at redshift z within a top-hat sphere. Furthermore, Sc{z) is the critical density contrast 
whose value predicted by the model for the collapse of a spherical top-hat fluctuation 
in a critical density universe, 6c = 1.69 independent of the redshift, has been confirmed 
by N-body simulations |^2||. In our analysis we used the expressions for 6c{z) provided 
in ref. for both low-density flat and open universes. We note, however, that at the 
redshift z=4.25, that we are considering, the resulting 6c value is always very close to 1.69. 

We note that the Press & Schechter approach [Q, on which eq. is based, holds 
only in the case of hierarchical clustering. In our case of WDM models, hierarchical 
clustering only takes place on mass scales larger than Mfs- On smaller scales, the lack 
of fluctuations causes the flattening of ctm- Therefore, by estimating ctm at arbitrarily 
small masses, one obtains the r.m.s. fluctuations at the free-streaming mass scale. In our 
approach, we will give up the dependence on mass scale M, which amounts to assume 
that DLASs are assumed to be hosted within protostructures of mass of about Mfs] 
protostructures of smaller mass, instead, are produced later by fragmentation of larger 
lumps. 

As for the observational value of Qhi, Storrie-Lombardi et al. provided for fio = 1, 
^Hi = (1-1 ± 0.2) X 10^^ h^^ a.t z = 4.25. In the light of all the above uncertainties in 
directly relating Qfji to ^Icoii, "we prefer to maintain a conservative approach here and to 
consider a model as ruled out if it predicts Qhi to be less than the observational la lower 
limit. At this level of comparison we do not consider as reliable to put constraints to 
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model producing too high a Qhi value. 

Furthermore, we should also rescale appropriately the value by Storrie-Lombardi to 
include the more general fio < 1 cases. Therefore, the limiting value that we consider is 



= 0.0009 h-'f{no, ^A, z = 4.25) , (23) 



where 



mo,n,,z) = \^^±j^^^^±^ ; n, = i-no. (24) 

3.2.2 The cluster abundance 

As for the cluster abundance, it has been recognized to be a sensitive constraint on the 
amplitude of the power spectrum [0. Based on the Press & Schechter approach |^ 



it is easy to recognize that the number density of clusters with mass exceeding a given 
value is exponentially sensitive to the r.m.s. fluctuation on the cluster mass scale. Fitting 
the local X-ray cluster temperature function with the Press-Schechter approach ||4^ led 



several authors to obtain rather stringent relationships between erg, the r.m.s. fluctuation 
value within a top-hat sphere of 8h~^Mpc radius, and Qq [Q. In the following we will 
resort to the constraint by Viana & Liddle, who provided the most conservative and, 
probably, realistic estimate of errors, mostly contributed by cosmic variance effects on the 
local cluster population: 

a(fio) = 0.36 + 0.31^0 - 0.28r]^ ; = 

a(fio) = 0.59 - 0.16^0 + 0.06fig ; = 1 - ^^o , (25) 

with uncertainties corresponding to the 95% confidence level. 

4 Discussion 

As for the purely WDM models, we plot in Figure 4 the constraints on the {Qq, h) plane, 
for g^, = 150, from DLAS and cluster abundance. Only scale-free primordial spectra (i.e.. 
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Hpr = 1) are considered here. Left and right panels correspond to the low-density flat 
(AWDM) and open (OWDM) cases, respectively. The solid line delimiting the coarsely 
shaded area indicates the limit for the region of the parameter space which is allowed 
by the observed Qhi in DLASs: model lying below such curves should be considered as 
ruled out, since they produce a too small Qhi value at z = 4.25. The cluster abundance 
constraint by eq.(^) is represented by the finely shaded region. The dashed curves connect 
models having the same age of the Universe: to = 9, 11, 13, 15 and 17 Gyrs from upper to 
lower curves. 

As a main result, we note that there is almost no overlapping between the regions 
allowed by the two observational constraints: for fixed values of the Hubble parame- 
ter, cluster abundance tends to select relatively smaller Qq in order to satisfy the low- 
normalization request of eg . (pSf ) . On the other hand, the DLAS constraint favour higher 
density parameters, which has the effect of both decreasing the free-streaming scale and 
to increase the small-scale power even in the absence of any free-streaming. Judging from 
this plot, one would conclude that the whole class of gravitino-dominated WDM models 
would be ruled out by combining constraints on the cluster and on the galaxy mass scale. 
A residual possibility seems to exist to reach a concordance for fio~ 0.4 (fio~ 0.5) and a 
high Hubble parameter, 1 {h^ 0.9) for OWDM (AWDM) models. However, two main 
problems arise in this case: (a) all the current determinations of the Hubble constant 
indicates 0.5 < h < 0.8 (b) the resulting age of the Universe would be definitely too 
small, especially for OWDM models, even on the light of the new recalibration of globular 
cluster ages, based on the recent data from the Hypparcos satellite |^ . 

We also checked the possibility of considering non-scale-free primordial spectra {upr ^ 
1), although results are not explicitly presented here. We verified that assuming either 
blue {n.p,. > 1) or red [upr < 1) spectra does not improve the situation. In the first case, 
power is added on small scales, with the result that smaller Qq are allowed by DLASs. 
However, the price to be paid is a rapid increase of the cluster abundance, that also pushes 
toward smaller Qq the finely shaded area. As for red spectra, the opposite situation occurs: 
the reduction of small-scale power leads both constraints to favour relatively larger Qq 
values, with no overlapping with the two allowed regions of the [Qq, h) plane ever attained. 

As a matter of fact, the situation becomes even worse when considering f2o = 1 
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WHDM models. Results for this class of models are reported in Figure 5 on the {Q,y, Upr) 
plane. Left and right panels are for h = 0.5 and 0.6, respectively; smaller h values are 
disfavoured by Hq determinations, while larger values are constrained by the age of the 
Universe. In both cases the regions allowed by DLAS and cluster abundance are largely 
disjoined, especially as higher Q^, are considered. Indeed, increasing the neutrino fraction 
has the effect of further reducing the power on small scales, thus further suppressing the 
high-redshift galaxy formation. 

Based on such results we should conclude that none of the variants of the WDM 
gravitino-dominated scenario is able to account at the same time for the relatively small 
abundance of clusters at low redshift and for the relatively high Qhi in collapsed structures 
at high redshift. It is worth reminding that this result has been obtained with the rather 
conservative choice of = 150. As we have shown in the previous section, more realistic 
value of g^, should be even smaller, thus putting WDM-dominated model in an even worse 
shape. 

Which are the consequences of such results on the low-energy SUSY breaking models 
that we described in Section 2? Of course, a first possibility is that gravitinos were so light 
as to be irrelevant from the point of view of cosmic structure formation. For instance, 
the current understanding of high-energy physics phenomenology would surely allow for 
rriQ ~ 1 eV. In this case, would be negligible. Of course, since G represents the LSP, 
the source for a cold DM component should be found in this case outside the spectrum 
of SUSY particles (e.g., axions).[] 

On the other hand, if a scenario with rriQ ~ 100 eV will turn out to be preferred, a 
non-negligible Qq can not be escaped. In this case, three possible alternative scenarios 
can be devised. The first one is to allow for cold + warm DM. However, since gravitinos 
have a much smaller free-streaming scale than neutrinos with rrii, ~ 5 eV, this scenario 
would suffer from the same pitfalls of the standard CDM one, unless one takes Qq < 1. 
The second possibility would be to have a substantially larger g^, so that gravitinos behave 
much like CDM. However, as we have seen in Section 2.2, it is not clear how a substantially 
larger g^, can be attained within plausible SUSY models. The third possibility would be 
to abandon the assumption of Gaussian fluctuations in favour of texture seeded galaxy 
^See however a recent proposal that a sneutrino in the messenger sector can be a CDM candidate Pq|. 
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formation [|^, which would ease the formation of high redshift objects. However, also 
this possibility has been recently shown to suffer from serious troubles in producing a 



viable power spectrum of density fluctuations |^0|], which make texture-based models as 
virtually ruled out. 

One may argue that the gravitino abundances will be diluted to a cosmologically 
negligible level by late-time entropy production. On the other hand, as the low value of 

suggests, the reheat temperature after the entropy production should be lower than 
the electro-weak scale to avoid the re-thermalization of the gravitinos, which severely 
constraints possible ways to generate the baryon asymmetry of the Universe. 



5 Conclusions 

In this paper we analysed the cosmological consequences of assuming the dark matter to 
be dominated by light gravitinos with mass in the range ~ 100 eV - 1 keV, as predicted 
by gauge mediated SUSY breaking (GMSB) models. We pointed out that gravitinos with 
such a mass behave like warm dark matter (WDM), since their free-streaming mass scale 
is comparable to the typical galaxy mass scale. 

After estimating the number of degrees of freedom of relativistic species at the grav- 
itino decoupling, g^, we resorted to a Boltzmann code to compute the appropriate WDM 
transfer functions. These are used as the starting point to compare gravitino-dominated 
model predictions to observational data about the abundance of HI within high-redshift 
damped Ly-a systems and about the abundance of local galaxy clusters. 

The main results of our analysis can be summarized as follows. 

(a) Low-density WDM models with both flat (ACDM) and open (OCDM) geometry 
can not satisfy the two observational constraints at the same time, unless a rather 
small ^0 value {"^ 0.4) and a rather large Hubble parameter 0.9) are assumed. 
However, such requests would conflict with measurements of the Hubble constant 
and with current constraints about the age of the Universe. 

(b) As for warm -|- hot (WHDM) models, we find that they have an even harder time. 
The combined free-streaming of both neutrinos and gravitinos generates a strong 
suppression of fluctuations at ~ 1 /i^^Mpc scale. This makes extremely difficult to 
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form high-redshift {z ~ 4) protogalactic objects if we require the model to match 
the I0W-2; cluster abundance. 

Based on such results we claim that no variant of a light gravitino DM dominated 
model is viable from the point of view of cosmic structure formation. Therefore, in the 
framework of GMSB models, this amounts to require the gravitino to be light enough 
{uiq^ 50 eV) so as to be cosmologically irrelevant (unless entropy production with a suf- 
ficiently low-reheat temperature dilutes the gravitino abundances). In this case, however, 
one would lose the LSP candidate for implementing a CDM-dominated scenario. 

As a concluding remark, we should point out that, from the point of view of the 
particle physics model building, we still lack an exhaustive construction of realistic GMSB 
schemes, in particular as far as the details of the messenger sector are concerned. In this 
respect we hope that our analysis may constitute a useful guideline for the intense work 
which is going on in the GMSB option. 
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Appendix 



In this appendix, we summarize the decay widths to gravitino which are needed in the 
calculation of in Section 2.2. We denote the gluino by g (with mass M3), the winos 
{U{2)l gauginos) by W^, (with mass M2), and the bino {U{1)y gaugino) by B (with 
mass Ml). We ignored possible mixing between the gauginos and higgsinos. 
The decay widths involving the gauginos are 

n9^S + G) - ^^^^ P6) 



V(W-^W- + G) ^ ^-p^^l-!5j (M,>™.) (27) 

T{W^^W^ + G) = ^^%rfl-^V (M,<mw) (28) 

r(B^, + G) = (29) 

ViW'^. + G) = (30) 



727r 




cos^ 9w 


Ml 


487r 




sin^ 9w 


M| 


487r 


m|M2, 


sin^ 


Mf 


487r 




sin^ 9w 


m| 


727r 


m|M|,; 


cos^ 9w 


M| 


487r 




cos^ 6*^^ 





i 

Mf 



r(s^z + G) = ^^z:2^(i-S) (Mi>m^) (31) 



r(Z^S + G') = ""7 ^2 7.2 I 1-^1 (Mi<mz) (32) 



V(W^Z + G) = ^";7.J^. (M.>m.) (33) 

r(^^^. + «, = (M,<™.). (34) 

where m^, m^^ are Z- and VF-gauge boson masses, respectively, and 9w represents the 
electro- weak mixing angle. 

The decay width of a slepton with mass mi to gravitino is given as 

ra^i + G)^- (35) 

^ ^ 487r m|M|; ^ ^ 

A similar expression is obtained for the decay width of a squark. 
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Figure captions 



Figure 1. The shape of the trasfer functions for WDM gravitino models. Left panel: the 
effect of varying g^, for i^o = 1 and h = 0.5; solid, dotted and dashed curves correspond to 
the CDM case, to g^^ = 200 and g^ = 100, respectively. Right panel: the effect of varying 
the Friedmann background; heavy and light curves corespond to the CDM and WDM 
with g^ = 200 cases, respectively. 

Figure 2. The shape of the transfer function for the warm + hot DM models. Left panel: 
the effect of varying g^ at a fixed value of ^2,^ = 0.25. Right panel: the effect of varying 
at a fixed g^ = 200. 

Figure 3. The mass-scale dependence of the r.m.s. density fluctuations within a top-hat 
sphere. Left and right panels are for the same models as reported in Figure 1. Heavy and 
light curves are for WDM and CDM cases. As for the WDM curves, the value of M at 
which they become lighter corresponds to the value of the free-streaming mass, defined 
according to eg. (pO]) . 

Figure 4. Observational constraints for COBE-normalized WDM models, with g^ = 150, 
on the (fio; h) parameter space, for fiat low-density (AWDM) and open (OWDM) models. 
The finely shaded area corresponds to the 95% c.l. region allowed by the cluster abun- 
dance, as estimated by Viana & Liddle |^ (see text). The heavy solid curve delimiting 
the coarsely shaded area represents the limit of the region allowed by the data about the 
Qhi in DLAS at z = 4.25, as given by Storrie-Lombardi et al. (see text); models 
lying below such curves are excluded. Horizontal dashed curves connect models having 
the same age of the Universe: to = 9, 11, 13, 15, 17 Gyrs from upper to lower curves. 
Figuree 5. Observational constraints for COBE-normalized WHDM models, with g^ = 
150, on the {Q^, np^) plane, for h = 0.5 (left panel) and h = 0.6 (right panels). A vanishing 
tensor mode contribution to CMB temperature anisotropics is assumed for Upr < 1 models. 
The two shaded areas have the same meaning as in Figure 4. 
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